Tunneling through a parabohc barrier viewed from Wigner phase space 
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We analyze the tunneling of a particle through a repulsive potential resulting from an inverted 
harmonic oscillator in the quantum mechanical phase space described by the Wigner function. In 
particular, we solve the partial differential equations in phase space determining the Wigner function 
of an energy eigenstate of the inverted oscillator. The reflection or transmission coefficients R or T 
are then given by the total weight of all classical phase space trajectories corresponding to energies 
below, or above the top of the barrier given by the Wigner function. 
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Tunneling [1] of a particle through a barrier is one of 
the striking phenomena of quantum mechanics [2]. In 
the special case of a repulsive quadratic potential, cor- 
responding for example to an inverted harmonic oscilla- 
tor [3] shown in Fig. 1(a), the transmission coefficient T 
takes the form [4] 
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depicted in Fig. 1(b). Here e = E/{hil) is the scaled 
energy which is the ratio of the eigenvalue E and the 
natural energy parameter hfl, where is the steepness of 
the quadratic barrier and h denotes the Planck constant 
divided by 2tt. 

The expression Eq. (1) has played a crucial role in the 
context of nuclear fission [5] . It usually emerges [5] from 
a semiclassical analysis [6, 7] of the Schrodinger equation 
of the inverted harmonic oscillator [3]. However, in the 
present Brief Report we rederive Eq. (1) from quantum 
phase space using the Wigner distribution function [8]. 
In particular, we show that Eq. (1) corresponds to the 
quantum mechanical weight of all classical trajectories [9] 
that have sufficient energy to go above the barrier. 

This result is counterintuitive since in the standard for- 
mulation [2] of quantum mechanics a la Heisenberg and 
Schrodinger an energy eigenstate does not contain ener- 
gies other than the eigenvalue. In contrast, the Wigner 
function [8] of such a state relies on the trajectories of all 
energies, however with positive or negative weights. 

We emphasize that the Wigner function of tunneling in 
the inverted harmonic oscillator has also been analyzed 
in Ref. [10]. The authors of this paper first derive the 
quadrature representation of the energy eigenfunctions 
and then perform the integral in the definition of the 
Wigner function. In contrast, we start from the two par- 
tial differential equations [11, 12] determining the Wigner 
function from phase space. Therefore, we find the Wigner 
function without ever going through the wave function. 
This approach is not only direct but also yields immedi- 




FIG. 1. Tunneling coefficient T of an energy eigenstate of 
eigenvalue E through a parabolic barrier (a) in its dependence 
on E (b) explained in terms of classical phase-space trajec- 
tories (c) subjected to the boundary conditions of a particle 
coming from the left (d). For three different energies - below, 
at the top of, and above the barrier - we depict the classical 
phase-space trajectories (c) which are either being reflected 
from, stopping at the top of, or going above the potential 
hill, respectively. The crossed line represents the separatrix 
in phase space separating the trajectories coming from the 
left and from the right. Hence, under normal scattering situ- 
ations only half of phase space is accessible depicted in (d) for 
a particle approaching from the left. The quantum mechani- 
cal transmission curve (b) is due to the quantum mechanical 
weight of all classical trajectories going above the barrier pro- 
vided by the Wigner function. 



ately the proposed interpretation of the tunneling coeffi- 
cient. 
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We study the tunneling of a particle of mass M through 
a quadratic barrier of steepness f2 expressed by the 
Hamiltonian 



2M 2 



(2) 



Here x and p denote the position and the coordinate of 
the particle. 

For this purpose we consider the Wigner function [8] 

W{x,p) dy e-^Py'"^ r (x - I) (x + I) 

(3) 

of an energy eigcnstate \E) oi H with wave function 
V'-E = V'-e(^)- However, instead of solving first the time 
independent Schrodinger equation HipE = Ei/je for ipE 
and then performing the integration in Eq. (3) pursued 
in Ref. [10], we analyze the partial differential equa- 
tions [11, 12] 
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and 
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for the Wigner function in phase space. We emphasize 
that Eqs. (4) and (5) are exact for the inverted harmonic 
oscillator. 

The classical Liouville equation (4) implies that We is 
constant along the classical phase space trajectories of a 
fixed energy H given by Eq. (2) and shown in Fig. 1(c), 
that is 



We{x,p) = WE/(hn) 



H{x,p) 
Ml 
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Next we take into account the boundary conditions 
associated with a scattering process. Two distinct possi- 
bilities offer themselves: (i) the particle approaches the 
barrier from the left, or (ii) it impinges from the right. 

The two cases manifest themselves in different classical 
phase space trajectories. Whereas the situation (i) is 
described by the trajectories in the domain above the 
separatrix 



p = MVtx, 



(7) 



depicted in Fig. 1(d), the case (ii) covers the area below 
it. 

Hence, for a particle coming from the left, the Wigner 
function w'j^^ of an energy eigenstate reads 

W^^\x,p) = We/^m,) {^^^) - ^^"^)' (8a) 



where O denotes the Heaviside step function. Hence, only 
the classical trajectories above the separatrix contribute 
to the Wigner function as shown in Fig. 1(d). 

Likewise, for a particle approaching from the right we 
find 



W'^>{x,p) - WE/(m) [ ^^-^ ] OiAinx - p). (8b) 



With the help of the familiar identity 
x5{x) = 



(9) 



for the Dirac delta function it is easy to verify that 
both expressions satisfy the Liouville equation (4) as 
long as the function is differentiable. The form 
of We — Weiif) corresponding to the scaled eigenvalue 
e = E / {Ml) in its dependence on the dimensional energy 
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of a classical trajectory is then determined by the 
Schrodinger equation (5) in phase space. Indeed, when 
we substitute the ansatz Eq. (8) into Eq. (5) we arrive at 
the ordinary differential equation 



V- 
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4(e - r])W, = 0. 
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Again we have made use of Eq. (9). It is remarkable that 
Eq. (11) is independent of the Heaviside step function. 
In order to solve Eq. (11) we make a Fourier ansatz 



(12) 



where the limits ti and T2 of the integration will be de- 
termined in a way as to simplify the differential equation 
for ~ w^{t) resulting from (11). Indeed, with the 
integral relation 



/ dr«;,(r)^il(e"'-) 
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and integration by parts we establish the identity 

vmiv) = T he(T2)e"'^= - u;e(Ti)e"'^^] 
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Similarly, we obtain 
d'^Ws 
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and 



dWe 
dr/ 



dr IT u'e(r) e"''^. 



2 „i?7Ti] 



(15) 



(16) 



which when substituted into (11) yields the ordinary dif- 
ferential equation 

i(4-T2)^-(ir + 4e)w, = (17) 
dr 

of first order for ~ w^{t). Here wc have made the 
choice Ti = — 2 and T2 = 2 which eliminates the boundary 
terms in Eqs. (14) and (15). 

By direct differentiation we can verify that 



Weir) = A{4-T^y 



' exp 
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is a solution of Eq. (17) and the constant 

TT 



(18) 
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of integration follows from the normalization condition 

f + 00 



1 
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which is a consequence of the Fourier ansatz Eq. (12). 

As a result, the crucial part of the Wigner function 
of an energy eigenstate of the inverted harmonic oscillator 
reads 



Ws(^) = 



1 



+2 exp 
dr 

-2 



-i£lnf|±^ 



(21) 



Since our ultimate goal is to calculate the transmission 
and reflection coefficients we do not discuss this expres- 
sion in more detail [13] but emphasize that is real 
and satisfies the symmetry relation 



(22) 



In Fig. 2 we depict the Wigner function of the en- 
ergy eigenstate of an inverted harmonic oscillator below 
the top of the barrier subjected to the initial condition 
of the particle approaching from the left represented by 
Fig. 1(d). In order to bring out the characteristic features 
we have rotated the phase space by 90° . The wave fronts 
in the foreground correspond to the phase-space trajec- 
tories of particles coming in and being reflected from the 
barrier. The waves on the upper left of the figure with 
much reduced amplitude represent particles going above 
the barrier. In the flat part above the separatrix extend- 
ing to the right upper corner the Wigner function van- 
ishes as to obey the boundary condition corresponding 
to the gray shaded area in Fig. 1(d). 

We now turn to the evaluation of the transmission 
probability T which is given by the integral 

POO 

T= / dT^WM (23) 
Jo 

of the Wigner function over positive values of rj only. 




FIG. 2. The Wigner function W^?' of an energy eigenstate of 
the inverted harmonic oscillator for the dimensionless energy 
e = —0.4 that is below the top of the barrier and correspond- 
ing to a particle coming from the left. Due to this boundary 
condition the Wigner function vanishes in the phase-space do- 
main above the separatrix (right upper part). It is constant 
along the classical trajectories with a dominant maximum at 
the trajectory corresponding to the energy eigenvalue (front). 
However, the Wigner function also reaches into the domain 
of trajectories traversing the barrier (left). The total weight 
underneath the Wigner function in this realm represents the 
tunneling coefficient. 



With the help of the identity 



dr/e'''^ = TrS{T)+iV 



we find from Eq. (21) the expression 



1 i 

2 47r 



— oo 



sinh(^/2) 



(24) 



(25) 



where we have introduced the new integration variable 
^ = In [(2 + t)/{2 — t)] and V denotes the Cauchy prin- 
cipal part. 

Due to the antisymmetry of sinh(^/2) only the imagi- 
nary part survives, that is 



1 



sin(e^) 



TT 7o sinh(^/2) 



and with the integral relation [14] 



sin(e^) n , / ett 
^^ihH^-2^^^^H^ 



we finally arrive at 



T = - [1 + tanh(£7r) 
in complete agreement with Eq. (1). 



(26) 
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FIG. 3. Reflection and transmission coefficients R = R{e) 
(left) or T = T(e) (riglit) in tlieir dependence on tlie energy 
eigenvalue e originating from the total Wigner function weight 
of all classical trajectories with energies rj below or above 
the top of the barrier, respectively. The symmetry relation 
Eq. (22) giving rise to i? = 1 — T is apparent. 



We conclude by noting that the reflection coefficient R 
follows from the identity 

R + T^l, (29) 

together with the normalization condition Eq. (20) and 
the definition Eq. (23) of T, as 

R= f dr^Weiv)- (30) 

J —OQ 

Therefore, it is the total quantum mechanical weight 
of the classical trajectories that are reflected, that is all 
energies below the maximum of the barrier. 

In summary, we have rederived the familiar reflection 
and the tunneling coefficients R and T for an inverted 
harmonic oscillator using the corresponding Wigner func- 
tion. This approach shows that R and T represent the 
quantum mechanical weights given by the Wigner func- 
tion of all classical trajectories that are reflected from or 



transverse the barrier as indicated in Fig. 3. The weight 
of a given trajectory follows from an ordinary differential 
equation of second order which we have solved using a 
Fourier representation of the Wigner function. 

Here we have not analyzed in detail the structure of 
the resulting differential equation of first order given by 
Eq. (17) in the complex plane. It suffices to say, that 
the origin of Eq. (1) is deeply rooted in the two poles at 
T = — 2 and r = +2 giving rise to a logarithmic phase 
singularity [15] contained according to Eq. (18) in the 
kernel of the Fourier representation of Wg . 

Such singularities also appear [3, 10] in the quadrature 
representation of the energy wave function and manifest 
themselves in the Unruh effect [16], the Hawking radi- 
ation [17] or optical analogues [18] of event horizons of 
black holes. To identify in the complex plane the crucial 
contributions to the integral Eq. (12), or to elaborate 
on the importance of the logarithmic singularity, and to 
compare and contrast the similarities and differences be- 
tween tunneling and particle creation goes beyond the 
scope of this Brief Report and will be the topic of a fu- 
ture publication. 
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